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Abstract 

We investigate the behaviour of a relativistic spherically symmetric radiative 
star with an accelerating, expanding and shearing interior matter distribution in 
the presence of anisotropic pressures. The junction condition can be written in 
standard form in three cases: linear, Bernoulli and Riccati equations. We can 
integrate the boundary condition in each case and three classes of new solutions 
are generated. For particular choices of the metric we investigate the physical 
properties and consider the limiting behaviour for large values of time. The 
causal temperature can also be found explicitly. 

1 Introduction 



The problem of radiative gravitational collapse was first investigated by Oppenheimer 
and Snyder-. Their interior spacetime is represented by a Friedmann-like solution for 
an isotropic homogeneous universe, and the exterior spacetime is described by the ex- 
terior Schwarzschild metric. The process of gravitational collapse is highly dissipative. 
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Therefore heat flow in the interior of the star must be present, and taken into account so 
that the interior solution of the radiating star can match to the Vaidya^ exterior metric 
at the boundary. The investigation of the gravitational behaviour of a collapsing star 
depends on the determination of the junction conditions matching the interior metric 
with the exterior Vaidya metric across the boundary of the star. Santos^ formulated 
the junction conditions for a shear-free fluid distribution with isotropic pressures and 
made it possibile to complete the model. His treatment paved the way to investigate 
physical features such as surface luminosity, dynamical stability, relaxation effects and 
temperature profiles. Raychaudhuri-, showed that the slowest collapse arises in the 
case of shear-free fluid interiors. Kolassis et al— assumed geodesic fluid trajectories 
when generating an exact model. Their model was generalised to include several new 
classes of solution in geodesic motion by Thirukkanesh and Maharaj-. In the past 
many investigations in radiating collapse have focussed on shear-free spacetimes with 
isotropic pressures (see the treatments of Herrera et air-, Maharaj and Govender-, 
Herrera et al— and Misthry et al—). 

The next stage of development was to include shear in the model of a radiating star. 
Naidu et al— included anisotropic pressures in the presence of shear for the interior 
spacetime and found simple exact solutions for geodesic fluid trajectories. This toy 
model was generalised by Rajah and Maharaj— by demonstrating solutions to a Ric- 
cati boundary equation governing the gravitational behaviour. The general situation 
requires a model which is expanding, accelerating and shearing. Noguiera and Chan—, 
modelling shear viscosity and bulk viscosity, attempted such a study but found that 
they needed to utilise numerical techniques to make progress. Some recent progress 
has been made in finding exact models for Euclidean stars by Herrera and Santos— 
and Govender et al—. In Euclidean stars with shear the areal radius and proper ra- 
dius are equal throughout the evolution of the radiating star. Our objective here is to 
show that it is possible to solve the relevant equations, for the general case, exactly in 
a systematic fashion. Our approach is the first analytic treatment to consider exact 
models with all the kinematical quantities present. We believe that these solutions will 
be helpful in studying physical features of a relativistic star in an astrophysical setting. 

In this paper we attempt to perform a systematic treatment of the governing equa- 
tion at the boundary of the relativistic star with the interior consisting of a fluid which 
has nonzero acceleration, expansion and shear. The junction condition is a nonlinear 
partial differential equation containing all three metric functions of spherical symme- 
try. In Section 2, we derive the field equations and the junction conditions. In Section 
3, we give the boundary differential equation governing the gravitational behaviour 
of a radiating, shearing and accelerating sphere. In Section 4, three new classes of 
exact solutions to the boundary condition are found in closed form. In Section 5, we 
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briefly investigate the physical features of the model generated and present the explicit 
form of the causal temperature for a particular choice of the metric functions. Some 
concluding remarks are made in Section 6. 



2 The model 

The most general form for the interior space time of a spherically symmetric collapsing 
star, which is expanding, accelerating and shearing, is given by the line metric 

ds^ = -A^de + B^dr^ + Y\de^ + sin^ Odcj)''), (1) 

where A, B and Y are in general functions of both the temporal coordinate t and the 
radial coordinate r. The existence of a fluid 4- velocity vector u enables us to introduce 
the kinematical quantities 

= ^t";6^^^ 6 = W^a, (Tab = ha^h^c-d), (2) 

where hat = gab+UaUb {habU^ = 0) is the symmetric projection tensor. The acceleration 
vector if' {ifUa = 0) represents the acceleration of the fluid particles relative to the 
congruences of u; the expansion scalar B measures the rate of increase of a volume of 
fluid element; the shear aab {cTabU^ = = cr"a) represents the tendency of a sphere to 
distort to an ellipsoid. For the comoving fluid 4-veloctiy = and the line element 
([T]), the acceleration vector -u", the expansion scalar and the magnitude of the shear 
scalar a are given by 

/A' \ 

(3a) 
(3b) 

(3c) 

where primes and dots on the metric functions denote differentiation with respect to r 
and t respectively. The energy momentum tensor for the interior matter distribution 
has the form 

Tab = (P + p)UaUb + PQab + QaUb + QbUa + T^ab, (4) 

where p is the density of the fluid, p is the isotropic pressure, is the heat flux vector 
and Tiab is the stress tensor. The stress tensor can be expressed as 

TTab = {Pr " Pt) (riaUb - ^^ab j , (5) 
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where pr is the radial pressure, pt is the tangential pressure and n is a unit radial vector 
given by = -^6f. The isotropic pressure 

P=^(Pr + 2pt) (6) 

relates the radial pressure and the tangential pressure. 

For the line element ([T]) and matter distribution (jlj) the coupled Einstein field 
equations become 

_ 2 BY 1 1 Y^ 
^ ~ A^BY^Y^^A^Y^ 

1 / Y" Y'^ B'Y'\ 

1 

Pr = -r. 



(7b) 



Pt 



where the heat fiux g° = (0, g, 0, 0) has only the nonvanishing radial component. A 
comprehensive treatment of the effects of anisotropy with heat flow in general relativity 
was carried out by Herrera et al-^; the flrst study with anisotropy appears to be in 
the treatment of Lemaitre~. The system of equations f l7a|) -f l7d|) governs the general 
model when describing matter distributions with anisotropic pressures in the presence 
of heat flux for a spherically symmetric relativistic stellar object. For this model f l7a|) - 
fl7dp describes the nonlinear gravitational interaction for a shearing matter distribution 
which is expanding and accelerating. From ( I7a|) -( l7d|) . we observe that if forms for the 
gravitational potentials A, B and Y are known, then the expressions for the matter 
variables p,pr,pt and q follow immediately. When the radial and tangential pressures 
are identical then pr = p± which generates an additional nonlinear partial differential 
equation called the condition of pressure isotropy. 

The Vaidya exterior spacetime- of a radiating star is given by 

ds^ = - (l - dv^ - 2dvdR + R\de^ + sin^ ~ 
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where m{v) denotes the mass of the fluid as measured by an observer at infinity. 
The hne element ([H]) represents coherent null radiation. The flow of the radiation is 
restricted to the radial direction relative to the hypersurface S, which represents the 
boundary of the star. The matching of the metric potentials and extrinsic curvature for 
the interior spacetime ([1]) and the exterior spacetime ([8]) produces junction conditions 
on the hypersurface S. These can be written as 



A{Rs,t)dt 
Y{Rj,,t) 

(Pr)s 




h 2—^ 

-Re dv 



dv. 



(9a) 
(9b) 

(9c) 

(9d) 



The junctions conditions fl^?a|) -f l9d|) were first derived by Santos'^ for a shear-free radi- 
ating relativistic star. It is important to note the nonvanishing of the radial pressure 
at the boundary S. Thus there is an additional differential equation (19d|) which has 
to be satisfied together with the system of Einstein field equations f l7a|) - fl7d|) . Junc- 
tion conditions similar to (19dp are important in describing phenomena which arise in 
astrophysics. Di Prisco et al— generated junction conditions relevant to spherical col- 
lapse with dissipation, in the presence of shear, for nonadiabatic charged fluids. Causal 
thermodynamics, in the context of the Israel- Stewart theory, was utilised by Herrera et 
al— to study viscous dissipative gravitational collapse in both the streaming out and 
diffusion approximations. 



3 The boundary condition 

Substituting ( l7b|) and (l7d|) in (l9d|) we obtain the boundary condition which has to be 
satisfied at the stellar surface. 

(A A'\ A 
- + —]YY + 2—YY' 
A B J B 

-2 A (^A' + b) YY' - l^y' + ^2 = 0. (10) 



Equation f lTOj) is the governing equation that determines the gravitational behaviour 
of the radiating anisotropic star with nonzero shear, acceleration and expansion. It 
is clear that flTUl) is highly nonlinear; it is difficult to solve without making certain 
simplifying assumptions. Some exact solutions to (|T0|) were found by Naidu et al— 
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and Rajah and Maharaj^^ for particles in geodesic motion (-it" = 0) but expansion 
7^ and shear cr 7^ 0. Chan^° considered the general case with 7^ 0, 7^ and 
shear a 7^ but no exact solutions were found. Instead the boundary condition was 
analysed numerically to study the physical features of the model, producing a final 
state where the star has radiated away mass during collapse. By assuming a relation 
between the metric functions B and Y for Euclidean stars Govender et al— found 
particular models with shear. 

Our intention is to solve f lTOj) exactly without restricting the functions. For conve- 
nience we rewrite (fTOj) in the following form 



B- 



Y 



y2 



AY' 2AYY' 



AY A 
~^Y' ^ 2YY' 



B' 



Y' 
Y' 



A' Y 
'AY' 



B + 



A' + 



AY' 
~TY 



0. (Ill 



In general f lTT]) is a Riccati equation in the gravitational potential B. This Riccati 
equation can be solved in special cases. 



4 Exact solutions 

The complexity and nonlinearity in f llip makes it difficult to solve in general. However 
particular exact solutions can be found if we view f lTT]) as a first order differential 
equation in the variable B and place restrictions on the bracketed expressions. We 
demonstrate this in the following three cases. 



4.1 Linear equation 

Note that equation ( JTTl) becomes a linear equation if we set 



Y Y^ 

+ 



A Y 



A 



AY' ' 2AYY' A^Y' 
This equation can be written as 



A- 



Y Y 
Y^TY 



A 



2YY' 



A^ 
2YY'' 



0. 



(12) 



(13) 



which is a Bernoulli equation in the variable A. Even though Y is an arbitrary function 
this equation can be integrated in general, and we have 

yy2 



A' 



(14) 



h{r)~Y' 

where h{r) is a function of integration. With the result (fl^ . we find that ( ITT]) becomes 

AY'' 



B 



Y' A' Y 

Y' " ~aY' 



B + 



A' + 



2Y 



0, 



(15) 



6 



which is hnear in B. 

The bracketed expressions in f[T^ contain the functions A, Y and their derivatives. 
In spite of this difficulty it is possible to solve (ITB]1 and obtain B in general. Therefore 
the solution for the junction condition (fTTj) can be given by 



A 



B 



Y 




(16a) 



dt X 



A 



77 + ^)exp 



A'Y 
AY' 



-dt 



dt 



(16b) 
(16c) 



where k{r) is a function of integration. We believe that fll6al) - fll6cp is a new solution to 
the boundary condition (fTOj) . Note that the gravitational potential Y{t,r) is an arbi- 
trary function in this class of solution. Once Y is specified then an explicit form for A 
is generated from (fT^ and the integrals in (I16ap - fll6cp can be evaluated. Consequently 
explicit forms for the metric functions A, B and Y can be found. The choice for Y 
should be made to provide a physically reasonable model. 



4.2 Bernoulli equation 

Observe that equation (fTT!) reduces to a Bernoulli equation if we set 

AY' 



Integrating this equation we get 



A' + 



Y 



2Y 

Cijt) 
A^ 



0. 



(17) 



where Ci{t) is a function of integration. Substituting f lTSj) into f lTTj) we obtain 



B 



3Ci_ A A' 
2C^~ A^A' 



B 



7 ACi 



^2 



+ 



A 



A^ 



2AA'Ci A^A' 2CiA' AA' ACfA' ACfA' 



B' 



(19) 



which is a Bernoulli equation in the variable B. 

The coefficients in ( fT9l) contain the functions A, Ci and their derivatives; however 
it can be integrated in general. On integrating ( |T9|) we can write 

A'C^^^ 



B 



A^[Jldt + gir)y 



(20) 
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where g{r) is a function of integration and for convenience we have defined 



7 Cy^ACi A'C^ 

i _|_ ^ i 

2 A^ 



2r^3/2 



CiCl^^ ACt 



+ 



3/2 



+ 



+ 



A' 



1/2' 



Therefore the functions 



A = A{t,r 
B 



A'Cf^^ 



A^[JIdt + gir)y 



A^' 



(21) 

(22a) 
(22b) 

(22c) 



satisfy the junction condition (fTTI) . The model (l22l) is an exact solution to the boundary 
condition flTUl) . Note that the gravitational potential A{t,r) is an arbitrary function in 
this class of solution. Once A is specified, together with the integration constants Ci, 
then an explicit form for / can be determined. Then the metric functions A, B and Y 
can be expressed in closed form in terms of elementary or special functions. The choice 
for A should be made on physical grounds. 



4.3 Inhomogeneous Riccati equation 

Note that equation (fTTI) has the form of an inhomogeneous Riccati equation if we set 



Y' ^ _ Q 

Y' ~ IY' ~ 



Integrating this equation we get 



A = Ya{t), 

where a{t) is a function of integration. In this case f llip becomes 



B 



Y{l + a^) 
2aYY' 



a 



a^Y' 



B' 



Y'a + 



YY'cL 
2Y 



(23) 
(24) 

(25) 



This is an inhomogenous Riccati equation which is difficult to analyse in general. How- 
ever we shall show that it is possible to integrate this equation by placing restrictions 
on the functions a and Y . 

If we take a to be a real constant and Y to be the separable function 



Y{t,r) = K{r)C{t), 



(26) 



where K{r) and C{t) are arbitrary functions of r and t respectively, then equation fl25l) 
becomes 



B 
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The Riccati equation (!27|l is not in standard form. Consequently we introduce the 
transformation 

B = wC (28) 



to obtain 



2aK' 



w = ^. (29) 



;i + a^)w^ - 2aK'w - 3a^K'^ 
The advantage of the form (l29l) is that it is a separable equation in the variables w 
and C. Equation ( l29l) can be integrated if the constant a is specified. To demonstrate 
a simple exact solution we take a = —2. Then (1291) can be written as 

(30) 



w 1 C 



{5w - 6K') {w + 2K') AK' C ' 

On integrating the above equation we obtain 

_ 2/f-|3C' + /(r)] 
f,C'-f(r) • 

where /(r) is a function of integration. The metric function B then follows since 
B = wC. 

Therefore we have generated a new solution to the inhomogenous Riccati equation 
(p5l) . The form of the solution is given by 

A = -2KC, (32a) 
_ 2K'C[3C^ + fir)] 

~ - m ' ^^^^^ 

Y = KC. (32c) 

The form of the solution fl32ap - fl32cp is particularly simple and does not involve further 
integration. The functions C, K and / are arbitrary; the physics of a specific model 
investigated will determine their explicit form. It is remarkable that the explicit so- 
lution f l32ap -( l32cll can be found for the boundary condition ( JTOi) in this case; Riccati 
equations are difficult to solve and only limited classes of solution are known to exist. 



5 Example 

The simple forms of the solutions found this paper make it possible to study the physical 
behaviour of the model. In this section we briefly consider the physical features of the 
solution generated in Section 4.3. For the gravitational potentials obtained in ([32a])- 
f l32cp . we take C{t) = t^, K{r) = r and /(r) = k, where is a real constant. For these 
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values the kinematical quantities become 



u 



e 



[5i-lf 



0,0 

4r[3t + l]Wkt 
[3 + 26t - 45P] 
2r[-l + 2t + 15t2](H) 

rTT' 
U 
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(33a) 
(33b) 
(33c) 



3r[l -2t- 15t2](H)- 

where we have set t = t^/k for convenience. From fl33ap -( |33cl) we observe that the 
acceleration m", the expansion B and the magnitude of the shear scalar are nonzero. 
These quantities remain finite in the interior apart from the stellar centre. Also note 
that in the limitine t oo the acceleration the shear scalar a — )■ 

and expansion — !■ 0. As the model evolves for large time the kinematical quantities 
grow progressively smaller. From the forms of ii"' and B given above we observe that 
the acceleration decreases more rapidly than the expansion for large time. 
The matter variables become 

[-3 - 43t + 15? + 95?] 



P 



Pr 



Pt 



2r2[3t + l]2[5t - r 
3[5t- 1]2 



kt 



Ar^Vkt V [3t + 1]^ 



4[-f -45t 



95? + 25? 



r^k[l -2t- 15t2]2 
[l + 25t- 165? + 75?] 



(34a) 
(34b) 
(34c) 
(34d) 

4r2[3t + l]3(Hy 

From fl34ap - fl34dp we observe that the energy density p, radial pressure pr, tangential 
pressure pt and heat flux q are continuous in the stellar interior, apart from the centre. 
At later times as t — )■ cxd we note that g — )■ so that the heat flux is radiated away during 
the process of gravitational collapse. It is interesting to see that the energy density 
p, radial pressure pr, the tangential pressure pt and the heat flux q are proportional 
to r~2^ and are decreasing functions as we approach the boundary of the star. The 
behaviour that p oc r"^ is of physical importance. It is interesting to observe that this 
property is also present in Newtonian isothermal spheres and relativistic isothermal 
cosmological models as pointed out by Saslaw et al—. 

A qualitative analysis of the matter variables, energy conditions, and stability is 
difficult to achieve for the interior matter distribution. However it is possible to gen- 
erate graphical plots which indicate physical viability. In Fig. [T]|3]we have plotted the 
energy density p, the radial pressure Pr and the tangential pressure pt- We observe 
that p > 0,pr > and pt > 0. In addition we have the behaviour p' < and p[. < so 
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that p and p are decreasing functions outwards from the centre to the stellar surface. 
For fixed values of the radial coordinate it is possible to plot the behaviour of 



Z = {p + p^f-Aq^ (35) 
Y = p-p,-2p4 + [(p + p,)2- V]f (36) 

Typical behaviour of these quantities are represented in Fig. [Hand Fig. |5] respectively. 
These graphs show that Z > and Y > 0. The behaviour exhibited in this physical 
analysis indicates that the weak, strong and dominant conditions are satisfied in interior 
points away from the centre. Also note from Fig. [6] that the speed of sound is less than 
the speed of light so that causality is not violated. 

Next we briefly consider the relativistic effect of the causal temperature in our 
model. The Maxwell-Cattaneo heat transport equation, in the absence of rotation and 
viscous stresses is given by 

rK% + qa = -K {hy,T + Tiia) , (37) 

where r is the relaxation time, k is the thermal conductivity, hat = Qah + UaU^ projects 
into the comoving rest space and T is the local causal temperature. Equation (l37j) re- 
duces to the acausal Fourier heat transport equation when r = 0. The causal transport 
equation (!37l) can be written as 

T(t,r) = — \- I \T{qB)B + AqB^] dr (38) 



kA _ 

for the metric ([1]). Martinez^, Govender et al^ and Di Frisco et al"^^ have shown that 
the relaxation time r has a major effect on the thermal evolution, particularly in the 
latter stages of collapse. Rajah and Maharaji2, and Naidu et al^ showed that in the 
presence of shear stress, the relaxation time decreases as the collapse proceeds and the 
central temperature increases. For our case, (138|) becomes 

r[-l + 15t-315P + 45P] 



T(t,r) 



Kr2[3t + l]2[5t- l]Vkt 
[l + 30t- 15?(H)3]ln[r] 



+ h{t), (39) 



Kr[-1 + 2t + 15t2](H)4 

where h{t) is a function of integration and we set r and k, as constant. When r = 0, we 
can regain the acausal (Eckart) temperature from (!39|) . It is possible to plot the causal 
and acausal temperatures against the radial coordinate. In Fig. [7] the temperature 
profiles are similar to the curves in Rajah and Maharaj-^. The temperature is a 
decreasing function from the centre to the boundary of the star in both the causal and 
acausal curves. The inclusion of particle acceleration in our model may contribute to 
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Figure 1: Density 
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Figure 2: Radial pressure 

the more rapid decrease of temperature from the core to the boundary of the star. This 
may be apphcable to phases of collapse where there is rapid expansion and cooling of 
the outer layers of the stellar fluid. As in the Rajah and Maharaj— model it is clear 
that the causal temperature is greater than the acausal temperature throughout the 
stellar interior. 

6 Discussion 

In summary, we considered the general case of a spherically symmetric radiative star 
undergoing gravitational collapse when the interior spacetime consists of an accelerat- 
ing, expanding and shearing matter distribution. The junction condition is rewritten 
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Figure 3: Tangential pressure 
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Figure 4: Z = {p + - Aq^ 
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Figure 6: Sound speed 
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so that it can be considered as a first order equation in the potential B{t,r). It is 
then possible to consider the junction condition as a standard differential equation: a 
linear equation, a Bernoulli equation and a Riccati equation. The linear and Bernoulli 
equations are solved in general. The Riccati equation can only be solved for a par- 
ticular value of the integration constant. Therefore three new classes of solutions to 
the boundary condition have been found. For a particular metric, corresponding to 
the inhomogenous Riccati equation, it is possible to obtain forms for the kinematical 
quantities and matter variables. It is then possible to indicate the behaviour of the 
model for large values of time. 
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